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Abstract- Eccentric cementing status always occurs in the period 
of well completion in complex formation. Under this condition, 
buckling failure of casing has greater impact on normal 
production. According to adjacent balance criterion and higher 
differential equation, a determination method of buckling load 
for eccentric cementing casing was established. This avoids 
complex modeling and discretization processes when determining 
the critical buckling load of casing, and geometrical and physical 
parameters can directly be put in the formula of critical buckling 
load. There exists 8.3% relative error, comparing critical load 
from the formula with the result by using finite element method. 
That means the results are accurate. This method has the 
advantage of direct calculation and is easy for engineering 
applications without complex modeling and discretization 
processes. 
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I INTRODUCTION 

Eccentric cementing means there appears eccentric 
annulus between the casing and borehole wall because the 
casing in the borehole is not centered. When cementing, the 
mud cannot be completely replaced by cement slurry. Part of 
cement sheath between wall and casing after the solidification 
of cement contacts with casing closely, and other part contacts 
with underground liquid. The state calls eccentric cementing 
status (Fig. 1), which often occurs in horizontal well section' 11 . 
It is the condition that does not meet the requirements of 
cementing. This condition often appears in well completion of 
complex formation . Obviously, eccentric cementing status 
lies in axis of casing. The casing and cement contact whit 
formation closely. On the ring of casing, only partially 
external circular surfaces of casing and cement contact with 
formation. Partially external circular surfaces of casing 
contacts with cement and formation, other meets with no 
borehole, annular space is full of formation liquid and mud. 
The load on casing is the hydrostatic pressure, which is the 
load of partially external circular surfaces by underground 
liquid (Fig. 2). The entire cement circle on the outer surface of 
casing could improve the collapsing strength' 31 . The paper [ ] ~ 
[8] studied the collapsing strength in common condition. In this 
paper, according to adjacent balance criterion, casing 
deformation of eccentric cementing status was studied and the 
analysis method for collapsing strength of eccentric cementing 
casing was confirmed by means of mathematical analysis 
solution . And the example was shown. 
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Fig. 1 Eccentric cementing chart 
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Fig. 2 Eccentric cementing mechanical model 



II. SIMPLIFICATION OF MECHANICAL MODEL FOR ECCENTRIC 
CEMENTING STATUS CASING 

Suppose there was considerable length for eccentric casing 
along axis (length of bare casing was beyond the range of 
constraint by cement), and it can be supposed as indefinite 
length. In order to simplify the issue, the physical dimensions 
of casing without cement were assumed: length as / along 
axis (considerable length), circumference size was expressed 
by centre angle of casing middle surface. That is 20 Q . 

Because the size of bare casing along axis was big enough, 
there was no influence between constraint along axis and 
circle. When studying the buckling deformation of eccentric 
cementing status casing, take unit length of casing (ring) for 
study. The ring could be replaced by the circle in which 
average diameter of casing lay. Mechanical model was showed 
as Fig. 3. 
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Fig. 3 Chart of mechanical model 

m. RELATION BETWEEN THE DEFORMATION OF RING 
AND THE LOAD 

As shown in Fig. 3b: unit length of the ring, average radius 
(the radius of middle surface) a , wall thickness t . 

Since — « 1 , the casing was thin-walled. Dots (Fig. 4) were 
a 

expressed with polar coordinates (r, 6) . 

Without deformation at ring, the geometrical relationship 

was: 



(1) 



x = r cos 6 
y = r sin 9 
z = r - a 
(ds) 2 = {rdOf 
According to energy principle, under the load of 
hydrostatic pressure P , deformation appeared at the ring, and 
the equation between deformations of local fixed ring and load 
1101 was: 
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In which: 

V ^ W — Displacement component of middle surface of the 
ring along circumferential and radial, m; 

E — Modulus of elasticity of the ring, Pa; 
A — Sectional area of the ring, m2; 

/ — Moment of inertia of section relative to middle surface, 
m4. 

IV. BUCKLING ANALYSIS OF THE RING 

Under the axisymmetric load, for the entire load P , there 
were two types of shape for the stability of ring. One was 
circular stability balance; the other was noncircular stability 
balance. Obviously, from circular stability balance to 
noncircular stability balance, there must be a balance in the 
form of critical state, which the load was corresponding to was 

called critical load, written as p cr . 

Assume circumferential and radial displacement 
component of circular stability balance as ( V , W Q ), and 
circumferential and radial displacement component of 
noncircular stability balance as ( V + Av,W + Aw ). 
Circular configuration and noncircular configuration both met 
Equation (2). In which ( Av,Aw ) was infinite small 

incremental. For configuration of circular stability balance, 
there was: 



v =v =w =o 



(3) 




Fig. 4 Thin-walled dots 



Then, take ( V,W ) to configuration of noncircular 
stability balance, (V, W ) could be showed as following. 



\V = Av 

[W = W + Aw 



(4) 



Take Equation (3) and (4) to Equation (2), and ( Av, Aw ) 
was considered. Omit high order, we got: 
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£A a 2 [(Av)'+(Aw)]' 

+ E/[(Av)-(Aw)']" 
-£A,^oPv)-(Aw)'] 

pa 3 [(Av)-(Aw)'] = 
E/[(Av)-(Aw)]'" 

£V 2 [(Av)'+(Aw)]" 

■£4 aW r [(Av)-(A W )'] 
- pa 3 [(Av) + (Aw)] = 
Due to symmetry: 

.2 



Assume: 



(5) 



W = 



pa 



(6) 



Take Equation (6) to Equation (5), the stability equation of 
local fixed ring impacted by hydrostatic pressure was 
obtained. 

A a 2 [(Av)'+(Aw)]'+ 
/[(Av)-(Aw)']"=0 

£V 2 [(Av)'+(Aw)] 
-E/[(Av)-(Aw)']"' 
+ pa 3 [(Aw) + (Aw)"] = 

According to adjacent balance criterion 1101 , ( Av , Aw ) in 
Equation (7) was circumferential and radial displacement 
increment of middle surface, which was from circular stability 
to noncircular stability of local fixed ring. Clearly, ( Av , Aw ) 
was the function of centre angle 9 of ring. According to the 
feature of structure symmetry and load symmetry, and 
( Av , Aw ) was the function of centre angle 6 , the function 
was the general solution for Equation (7), and it must meet the 
following boundary conditions. 



Av 
Aw 



-*„ = (Av) 
-a. = (Aw)' 



= 



= 



(8) 



(Aw)'L =0 =0 



Suppose 



a l = EA a 2 ,a 2 =EI,a 3 = pa' 

[X = (Av)'+(Aw) 
[Y = (Av) - (Aw)' 

Equation (7) can be rewritten as: 

( ai X'+a 2 Y"=0 

\a x X -a 2 Y"'+a 3 [X -F] = 



(9) 



(10) 



To Equation (7) and (10), analytic solutions can be 
obtained with high differential equation method. 
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Equation (10) can be rewritten as first-order equations. 
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In which A was 
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(15) 
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Characteristic polynomial for Equation (14) 

\Xl + A\ = A 2 (A 2 + + = 

Characteristic value for Equation (16): 



(16) 



X x = 0(doub 1 e root) 



K = 



a 3 
a, 



(17) 



a 3 
a. 



So that three equations were acquired as following. 

(-A-A l l) 2 u=0 
< (-A-A 2 l)u=0 
(-A-l 3 l)u=0 



(18) 
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Take tl l — 2,n 2 — l,n 3 — 1 , and the basic solution for 
Equation (18) was: 



In which 
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(19) 



Suppose 77 = (^,^ 2 ,77 3 ,77 4 ) r =4+^,+^, , the 
relationship among GC,j3,/,T was 



a 2 a 3 ?] 3 - a 2 (a x + a 3 )rj x 
a x a 2 + a 2 a 3 + a 3 a x 



J3 = %~ 

_ rj 4 {a x a 2 +a 2 a 3 + a 3 a) + a 1 a 3 77 3 /l 2 - a x (a { +a 3 )?] l A 2 



r = 



2/l 2 {a x a 2 + a 2 a 3 + a 3 a x ) 



VA^ a i a 2 + a 2 a 3 + a 3 a ) )-a 1 a 3 rj 3 A 2 + a l (a 1 + a 3 )rj x X 1 
2/^ 2 (a Y a 2 + a 2 a 3 + a 3 a x ) 

Take Oi,f3,y,T to Equation (19), the expression of 
, <^ 2 , <^ 3 can be obtained. 

The form of solution for Equation (14) was 



(20) 



Assume fj = (l,0,0,0f . Take (0,1,0,0) T , 
(0,0,1,0) T and (0,0,0,1) T to (pit) , and four linearly 
independent solutions (p i (t) were got. Then from Equation 
(11) 

z4KYJjJ4^t)4t)4t)^(t))c 



(21) 



In whichC = (c 1 ,C 2 ,C 3 ,C 4 ) r , constant vector. So 



X=a l a 3 q +[a 2 <^ -a^ cos//] % 



Y=a l (a l -+a)tc l +c 2 + 



a 2 (a,-KyfH — —an// 



C3+(l+cos//)c 4 

(22) 



Equation (9) can be rewritten as 
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(23) 



Equation (23) was nonhomogeneous linear system of 
differential equations. The proper vector corresponding to the 

proper value ( = i , A 2 = — i ) of coefficient matrix 

'0 -1 A 



1 



That is: 



met the following equation. 

(I/-a)J=o 



(24) 
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(25) 



v 1 z 
Solution vector corresponding to it was 

ft - i 



(26) 



(27) 



and 



£1 

r 2 



- ; 
1 



(28) 



Sor ] =(/,l)' ,r 2 =(-/,l)' , and basic solutions for Equation 
(23) were 



<b(e)= 



K e e j 



(29) 



Inverse of Equation (29) was 



2 
1 



(30) 



The form of general solution for nonhomogeneous linear 
system of differential equations (Equation (23)) was 
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X(s) 
-Y(s) 



ds 



V v J J 
(31) 

What can be obtained according to the first and second 
term of Equation (8) was 



v Aw) ^\w(0-s) a^d-s) A-Y(s)J 



ds 
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From Equation (32) 



aa 

a 2 (a^)0 o A sm/y^ c 3 -(1-kes/^)c 4 



(34) 



Aw = f g [X (s) sin(0 -s)-Y(s) cos(0 -sj]ds 

The expression for function X(s)andY(s) can be seen 
in Equation (22). 

What can be obtained according to the first and second 
term of Equation (8) was 

4[-^sm^^(q^ocB^+(^^)^^-(^+^)^^]q 



(35) 



Third-order equations consist of four constant 
coefficients C 1 , C 2 , C 3 , C 4 : 



0=ajC^q -f{ a ^ ccs/Z^,) <^ -^sm/Zjjc, 



C,-(1-K3CS/^)C 4 



0=a,(c( 1 -(^)^q-q+ 

0=[-q 2 sin^+q(q -^axf^q -Hl-ax^ 

^ 2 2q (/*1) 2 2q (w-1) 4 

(36) 

Or be written as 



A n Cj + A 12 c 2 + A 13 c 3 + A 14 c 4 — 
AjjC, + A 22 c 2 + A 23 c 3 + A 24 c 4 = 

"^31^*1 ^32^*2 J ^33^ > 3 ^34^4 ~ ^ 



(37) 



In which 



4 ^flft,4 =Q4 -ofl,cEs/^,4=-^sin/^ 



4 ^(a,^^ =-1,4 ^a/^-to,)^ 4-^ sin/^,4 =-( 1-kes/^) 



a t q . 
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4 =-q 2 sin^-hq(q -ha^asf^ =l-ccs$ 

4=^^^(^«^^)^^^)-^ 



2^-1) 



q 2^fl) q 2T//-1) 
What was obtained was the proportional relationship among 
the four coefficients. Firstly, make C 4 = 1 , so Cj and C 3 were 
acquired from the first and third equation. Then take them to 

the second equation, we got C 2 . The expression of 

was acquired by taking them to Equation (32). 
On the basis of analysis, Av and Aw were the function 

of E , Aq , a , / , p and , it meant they were the function of 

Oj , a 2 , a 3 and in analyzing. Under the axial symmetry 

load, to make the ring be in critical state from circular stability 
balance to noncircular stability balance, there must have the 
following equations for all . 



|Av = 
Aw = 



(38) 



Meanwhile, critical load p cr of buckling failure of ring 
was gained. 

V. EXAMPLES 

Finite element model of eccentric cementing status was 
showed in Fig. 5. 




Fig. 5 Element analysis model of structure 

Basic parameter of casing: grade of steel was N80; 
external diameter (D) was 139.7 and wall thickness was 7.72; 
modulus of elasticity was 216GPa; poisson's ratio was 0.25. 
Stratum parameter: modulus of elasticity was 216GPa; 
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poisson's ratio was 0.3. angle of bare pipe O equaled to 
120° symmetrically. The load was external pressure P . 

Take the basic parameters above to the coefficients of 
Equation (7) and (9), we got: 

a = 69.85 xl(T 3 (m) ; A = 7.72xl(T 6 (m2) ; 
/=3.8342xl(T n (m4); Oj =EA a 2 =7909.8713 (Nm2) ; 
a 2 =EI = 0.8052(Nm2) ; a 3 = 0.06985/? (N/m2) 

After taking a i ,a 2 Q 3 and$ to Equation (36), when 

given P an initial value p Q , the first maximal value p cr of P 
was 155.1 MPa through trial. 

Under the condition above, the hydrostatic pressure 
P cr was 155.1 MPa when theoretical eccentric cementing 
appeared stability disruption. 

Finite element method: 

According to the physical model established with the 
parameter above, take plane 182 for discretization to the 
structure. Make horizontal displacement of left and right side 
limited, as well as vertical displacement of underside. Casing 
would be affected by external pressure p . 




(a)100M (b)200MPa (c)300MPa 

Fig. 6 Deflection of eccentric cementing casing 



There existed no relative displacement between casing and 
layer. Finite element model could be seen in Fig. 5. There 
were 1780 elements and 2008 nodes in total. 

In Fig. 6a to 6c, structural deforming graphs were showed 
under the condition of external pressure lOOMPa, 200MPa and 
300MPa. It can be seen from the graphs that both Fig. 6b and 
Fig. 6c were in unstable state. The relationship between radial 
displacement and external pressure on the peak of casing was 
given by Fig. 7, from which the break of displacement 
appeared when external pressure was 168MPa. That was to 
say, critical pressure was 168Mpa when it was unstable. 




Fig. 7 Radial displacement 
vs. pressure of casing at acme 

Compare finite element results with theoretical results, 
168-155.1 



error of calculation was: . 



155.1 



.3% 



Theoretical results 



were close to finite element results. And the accuracy of 
theoretical analysis was validated. 

VI. CONCLUSION 

The method, buckling analysis of eccentric cementing 
status casing, is an easy and valid way for the study of 
collapse resistance. It is not as complex as the physical finite 
element modeling process, but the stable load can be acquired 
by directly entering the relevant parameters. For the 
convenience of calculation, accuracy and convenience for 
engineering applications, it is easy for theory to be applied to 
engineering project. 



REFERENCES 



[1] 



Xu Huifeng. Drilling Technical Manual (3) (Cement) [M]. Petroleum 
industry Press, Beijing: 1990. 7. 

Wu Jiang, Hao Junfang. Permit eccentricity of well cementation casing 
[J]. Petroleum Drilling Techniques, 1989, 17(4): 42-47. 

Song Ming, Influencing Tendencies of Cement Sheath on Casing 
Retaining Capacity [J], Oil Drilling & Production Technology, 2002, 
24(4):7-9. 

Cui xiaobing, Zhang Hong, Li Zihou, An Wenhua, Song Zhi. Research 
on Reasonable Free Segment Length of Deep Oil Well Casing String [J], 
Acta Petrolei Sinica, 1997, 18(1): 91-95. 

Zhang Zhiquan, Wu Jiazhong, Luo Yue. Analysis and Countermeasures 
of Casing Deformation Causes in Gangxi Oil field [J], Journal of 
Jianghan Prtroleum Institute, 2002, 24(4): 83-84. 

Fang Jun. Analysis of Loading Propekty of Casing and Cement Sheath 
Under Nonuniform Geologic Stress [J], Journal of the University of 
Petroleum, 1995,19(6):52-57. 

Gao Xueshi, Fang Jun. Calculation Method for Collapse of Oil Casing 
Restrained by Concrete [J]. Journal of the University of Petroleum, 1994, 
18(3):39-46. 

Ma Xu, Lab Test of Destruction of Cement Sheath Outside Casing, 
China Petroleum Machinery, 2000, 28(4): 15- 18. 

Li Weiguo, Tong Dengke. Numerical Computing Methods [M]. 
University of Petroleum Press, Shandong: 262-291(2008.10 ). 

[10] Chen Tieyun, Shen Huizhong. The Bucking of Structure [M]. Shanghai 
Science and Technology Literature Press, Shanghai: 1 1 1-120(1993, 12). 



[21 



[31 



[4] 



[51 



[61 



[7] 



[81 



[91 



JPSR Vol. 1, Iss. 2, Jul. 2012 PP. 20-25 http://www.jpsr.org © American V-King Scientific Publishing 

-25- 



